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Данная работа посвящена описанию модифицированного алгоритма бинаризации цветных фотоснимков 
горных пород с целью выделения объектов интереса, являющихся различными по величине фракциями 
минералов, представленными на неоднородном фоне, для обеспечения в дальнейшем возможности 
автоматизации выполнения гранулометрического анализа их состава. Представлены результаты эксперимен- 
тальных исследований возможностей предложенного алгоритма, подтверждающие повышение достоверности 
бинаризации и выделения контуров фракций минералов по сравнению с исходным методом. 
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Введение 


Одним из основных этапов технологического процесса в горной промышленности 
является определение гранулометрического состава горных пород. Оно представляет 
собой расчет процентного весового содержания в породе различных по величине 
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фракций [1]. Для решения этой задачи традиционно применяется ситовый метод, 
основанный на просеивании породы через набор различных по величине ячеек [2]. 
Этот метод требует существенного применения ручного труда и временных затрат, 
что затрудняет его применение в ряде случаев, например, для проведения экспресс- 
анализа. Таким образом, актуальной является задача разработки технологий, позво- 
ляющих осуществлять определение гранулометрического состава горных пород в 
автоматизированном или автоматическом режимах с минимальными затратами ручного 
труда и обладающих достаточным уровнем точности. 

Одной из возможностей решения задачи проведения экспресс-контроля грануло- 
метрического состава является использование фотографических изображений, полу- 
ченных в процессе добычи или переработки [3] горных пород. При использовании 
такого подхода одним из важнейших этапов, определяющих точность конечного ре- 
зультата, является предварительная обработка исходного снимка, которая является 
нетривиальной задачей [4]. Этот этап включает в себя такие операции как фильтрация 
шумов, повышение контраста, выделение объектов интереса (контуров), бинаризацию и т.д., 
способ выполнения которых существенно зависит от качества и характеристик 
исходных изображений. Конечной целью обработки является получение данных в 
виде, пригодном для автоматического (автоматизированного) анализа гранулометри- 
ческого состава. 

Основными проблемами, усложняющими предварительную обработку изобра- 
жений, являются неоднородность фона и освещения при формировании снимка, 
наличие размытых и низкоконтрастных областей, «слипание» и частичное перекры- 
тие (наложение) объектами друг друга. 

Так, например, неоднородность фона и освещения затрудняет применение тра- 
диционных способов фильтрации (медианного фильтра или фильтра Винера [5]) и 
улучшение разделения объектов на основе применения гистограммных методов [5-7] 
или преобразования растяжения контраста [6]. Наличие низкоконтрастных и засве- 
ченных областей делает невозможным непосредственное выделение контуров объектов 
интереса на основе использования пороговых (метод Хуанга и Ванга [8]) или гра- 
диентных (оператор Собеля [5], [6]) методов. 

Постановка задачи. В работе [9] был предложен алгоритм бинаризации фотосни- 
мков горных пород позволяющий выделять контуры объектов интереса с целью прове- 
дения дальнейшего гранулометрического анализа. Его недостатком является выделение 
значительного количества неинформативных (ложных) контуров, отражающих струк- 
туру, а не границы объектов интереса, а также недостаточная чувствительность при 
разделении соприкасающихся объектов с одинаковыми яркостными свойствами, что 
обусловлено использованием только одного (красного) цветового канала. 

В данной работе предложен модифицированный алгоритм бинаризации цветных 
фотоснимков горных пород, который за счет учета всех цветовых каналов позволяет 
повысить достоверность их анализа. 


Решение задачи 


Предложенный алгоритм состоит из следующих шагов: 


1 
а) повышение яркости и контраста красного цветового канала /[, (выбран как 
наиболее информативный на основании экспериментальных исследований) после 
масштабирования его на отрезок [0,1]; 


6) удаление неоднородного фона (следует отметить, что его выбор оказывает 
существенное влияние на способ обработки изображения); 
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в) преобразование яркости полученного изображения с целью улучшения раз- 
деления объектов; 

г) выделение контуров; 

д) бинаризация полученного изображения; 

е) синтез изображения на основе индексов максимального значения цветовой 
составляющей каждого пикселя; 

ж) подавление малозначимых различий у сформированного изображения; 

3) выполнение шагов в) — г). Следует отметить, что вид преобразований, при- 
меняемых на этом этапе, несколько отличается от формул, используемых на шагах в)-гГ); 

и) шаг д) для изображения, полученного после выполнения предыдущего этапа; 

к) формирование итогового изображения на основании слияния двух бинарных 
снимков, полученных на шагах е) ии). 


1 
Повышение яркости и контраста для изображения 1 ‚ выполняется по формуле: 


(о НТ (0 


г 


для всех пикселей изображения, яркость которых удовлетворяют условию 


где | = [п + 0.5/2 ‚а 1 — среднее по изображению 1 

Для нейтрализации (удаления) фона осуществляется преобразование яркости 
каждого пикселя изображения /” по следующему алгоритму. 

1. Изображение [” разбивается на перекрывающиеся окна размерностью 3х3 
пикселя (во всех оконных преобразованиях предложенного алгоритма применяются 


окна такого размера), и для каждого пикселя \, текущего окна выполняется пре- 


и 
образование: 


ио=и +, (2) 


х,У х›У 2 


3 
где у? — его средняя яркость, в результате чего формируется изображение Г. 


2. Получение изображения /“ на основе степенных преобразований, обеспе- 
чивающих нелинейность преобразования яркостей: 


ное, © 


Этот шаг является одним из основных При подавлении (приведении к белому 


цвету) фона. 
3. Последовательное двукратное применение ранговой фильтрации (по макси- 
муму для области 3х3 пикселя) [10] с последующим выполнением 1-го шага, в ре- 


5 
зультате чего формируется изображение /`. 


4. Получение изображения /[ 6 по следующей формуле: 


о 4 


ху 
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5. Формирование изображения /[ й следующим образом: 


в Е (“ 


(5) 
6. На последнем шаге пиксели происходит последовательная замена всех зна- 
чений яркости, которые превышают М ‚› на 1. Это значение рассчитывается как 
среднее по яркостям изображения, удовлетворяющим условию: 
а (6) 
1 х,у : 


Этот процесс происходит до тех пор, пока выполняется условие М > 1+а-т)/З. 


На основе экспериментальных исследований рекомендуется использовать 11 = 0.9 
приравниваются 1. Затем применяется ранговая фильтрация (по максимуму) в ре- 


8 
зультате чего и формируется изображение /°. 
Преобразование яркости с целью улучшения разделения объектов выполняется 
с использованием следующих формул: 


в, = мы 50 п) 


НР в. 

а | (8) 

где /° — среднее значение яркости для пикселей изображения /“, чья интен- 

сивность меньше 1. После этого выполняется сглаживание на основе формулы (2), в 
результате чего и формируется изображение /Г'. 

Для выделения контуров использовался лапласиан гауссиана [6] со значением 

о =0,5 ‚ в результате расчета которого формировалась маска размером 3х3 пикселя. 


После применения полученной маски к снимку Г' формируется изображение /”, для 
которого вычисляется минимальное значение 1. Т.к. одним из недостатков приме- 
нения лапласиана гауссиана является формирование двойных контуров, то для всех 


> 12 
пикселей изображения /„,, значения яркости которых меньше 0, они обнуляются. Далее 


ши ? 


ВЫПОЛНЯЮТСЯ преобразования, усиливающие различия между цветом контура и фона: 


мы т 
(2 р пив _ 712 

Е = (1 ху та 9 
а > (9) 

13 

1-7 

Е: ху 

Ту — (из, у (10) 
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Производится бинаризация полученного изображения /“ с использованием 
условия 1. >71‘. тде Г” = то среднее по пикселям снимка / “ значения интен- 
сивности которых больше 0. Затем выполняется ранговая фильтрация порядка 7 для, 
в результате чего формируется первое бинарное изображение / ы я 


Для улучшения разделения объектов с различными яркостными и цветовыми 
характеристиками необходимо использовать информацию обо всех его цветовых ка- 
налах исходного изображения. В предложенном алгоритме это выполняется следую- 


щим образом: формируется изображение /', на основе индексов максимального значе- 
ния цветовой составляющей каждого пикселя и масштабируется на отрезок [0,1]. 

После этого осуществляется подавление малозначимых различий между значе- 
ниями яркости одного и того же пикселя у разных цветовых каналов. Это осущест- 
вляется путем присвоения 1 значениям яркости тех пикселей у изображения /[ а 
для которых разница их интенсивностей у любых двух выбранных цветовых каналов 
не превышает значения 4, / 4. При этом гы — среднее по разности между интенсив- 


ностями для пары выбранных цветовых каналов. Затем происходит медианная фильт- 
рация полученного снимка и дважды выполняется ранговая фильтрация (по максимуму). 
Полученное изображение подвергается преобразованию по формуле (2) в 


2 
результате чего формируется снимок /[;,/. Преобразование яркости изображения 


2 
[1 с целью улучшения разделения объектов осуществляется по формуле: 


(р 
( а } мы ( м }. $ 5. |] м, (11) 


после чего дважды выполняется ранговая фильтрация (по максимуму) для области 
3х3 пикселя. 
3 
Выделение контуров у изображения /[ а Также осуществляется на основе при- 
менения лапласиан гауссиана с вышеописанными параметрами. После применения 
полученной маски к снимку /[ . „; формируется изображение / - 


т › 


для которого вычисляет- 
ся минимальное значение /“., 
обнуляются. Далее осуществляется усиление различия между цветом контура и фона на 
основе преобразования: 


(5, = (№). (4) " а 


и формулы (10) в результате чего формируется изображение / о : 


И ДЛЯ пикселей, значения яркости которых меньше 0, ОНИ 


14 
Бинаризация этого изображения осуществляется также как и для снимка [”. В 
2 
результате этого формируется изображение /[,,. 


Итоговое изображение /“”" является результатом слияния снимков Г) и [., 


(на основе максимума яркости) и последующего применения ранговой фильтрации 
порядка 4. 

Экспериментальные результаты. Представленный модифицированный алгоритм 
был применен для обработки изображений, содержащих различные минералы горных 
пород на неоднородном фоне, примером которого может служить снимок, цветовые 
каналы которого представлены на рис. | а) — 1 в). 
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Рисунок 1 — Цветовые каналы исходного изображения: 
а) красный; 6) зеленый; в) синий 


На рис. 2 а) представлен результат повышения яркости и контраста красного 
цветового канала в результате чего, в частности, наблюдается и улучшение визуаль- 
ной различимости камня в левом верхнем углу. На рис. 2 6) приведены результаты 
удаления фона (замены его белым цветом), а на рис 2 в) — улучшения разделения объек- 
тов. Последний этап обеспечивает также уменьшение неинформативной детализации 
частей объектов, обусловленной особенностями их формы и структуры. На рис.2 г) 
отображен результат улучшения разделения объектов на изображении, которое было 
получено на основе индексов цветового канала с максимальной яркостью для каждого 
пикселя. Это изображение позволяет лучше разделить объекты различных цветов. 

На рис. 3 а) представлены результаты бинаризации красного цветового канала для 
исходного алгоритма. На рис 3 6) — 3 г) продемонстрировано применение модифициро- 


ванного алгоритма бинаризации. При этом на рис. 3 6) приведено изображение / ры (би- 


2 
наризация красного цветового канала), на рис. 3 в) — снимок /р.„, а на рис. 3 г) — конеч- 


ный результат. Очевидно, что использование информации о цветовых каналах позволило 
улучшить разделение объектов различных цветов. При сравнении результатов бинариза- 
ции красного цветового канала исходного (рис. 3 а) и конечного результата модифици- 
рованного (рис. 3 г) алгоритмов можно увидеть, что в последнем случае уменьшается 
количество выделенных ложных контуров, а также улучшается разделение объектов. 
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в) г) 
Рисунок 2 — Результаты обработки исходного снимка: а) повышение яркости и контраста для 
красного цветового канала; 6) удаление фона; в) улучшение разделения объектов 


и). 3 
(изображение /); г) улучшение разделения объектов (изображение /.„) 


Рисунок 3 — Результаты бинаризации: а) применение исходного алгоритма; 6) изображение 


1 2 
р 5в- изображение бя $5 Г - итоговое изображение 
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Выводы 


Предложенный в данной работе модифицированный алгоритм бинаризации цветных 
фотоснимков горных пород с целью выделения объектов интереса для проведения 
дальнейшего гранулометрического анализа позволяет повысить достоверность их 
выявления по сравнению с исходным методом. К недостаткам предложенного подхода 
следует отнести выявление, хоть и в меньшей степени, ложных контуров, обуслов- 
ленных структурой объектов, невозможность обработки засвеченных объектов (что 
выражается в интерпретации нескольких таких объектов в качестве единой струк- 
туры) и большую вычислительную сложность. 
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